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ABSTRACT 

We present a new sample of purely near infra-red selected i^vega <16.5 
[^AB <18.4] extremely red ((J — i^)vega > 2.5) quasar candidates at z ~2 from 
:::^900deg^ of data in the UKIDSS Large Area Survey (LAS). Five of these are spec- 
troscopically confirmed to be heavily reddened Type 1 AGN with broad emission 
lines bringing our total sample of reddened quasars from the UKIDSS-LAS to 12 at 
z=1.4-2.7. At these redshifts, Ha (6563A) is in the iC-band. However, the mean Ha 
equivalent width of the reddened quasars is only ten per cent larger than that of the 
optically selected population and cannot explain the extreme colours. Instead, dust 
extinction of Ay ~2-6 mags is required to reproduce the continuum colours of our 
sources. This is comparable to the dust extinctions seen in submillimetre galaxies at 
similar redshifts. We argue that the AGN are likely being observed in a relatively 
short-lived breakout phase when they are expelling gas and dust following a massive 
starburst, subsequently turning into UV-luminous quasars. Some of our quasars show 
direct evidence for strong outflows (v~800-1000 km/s) affecting the Ha line consistent 
with this scenario. We predict that a larger fraction of reddened quasar hosts are likely 
to be submillimetre bright compared to the UV-luminous quasar population. We use 
our sample to place new constraints on the fraction of obscured Type 1 AGN likely 
to be missed in optical surveys. Taken at face- value our findings suggest that the ob- 
scured fraction depends on quasar luminosity. The space density of obscured quasars 
is ~5x that inferred for UV-bright quasars from the SDSS luminosity function at 
Mi < —30 but seems to drop at lower luminosities even accounting for various sources 
of incompleteness in our sample. We find that at M^ ~ —28 for example, this fraction 
is unlikely to be larger than ~20 per cent although these fractions are highly uncertain 
at present due to the small size of our sample. A deeper if-band survey for highly 
obscured quasars is clearly needed to test this hypothesis fully and is now becoming 
possible with new sensitive all-sky infra-red surveys such as the VISTA Hemisphere 
Survey and the WISE All Sky Survey. 

Key words: galaxies: active, (galaxies:) quasars: emission lines 



1 INTRODUCTION 

The last few years have seen an explosion in multi- 
wavelength surveys covering very large areas of sky and 
extending to wavelengths from the X-ray through the 
optical/infra-red and into the sub-millimetre, millimetre and 
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radio. Such surveys are not only allowing us to gain a more 
holistic view of galaxy and active galactic nuclei (AGN) 
properties by sampling their entire spectral energy distri- 
bution (SED), but their unprecedented area also opens up 
discovery space for rar e and unu sual classes of objects e.g. 
the most distant (iMortlock et al.« 20 11) as well as the most 
luminous quasars (|lrwin et al .11 19981 ). Observational studies 
of distant quasars and galaxies have gained new impetus re- 
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cently with the advent of sensitive infra-red detectors which 
are allowing us to observe the SEDs of galaxies and quasars 
after they have been redshifted beyond optical wavelengths. 
In addition to sources which appear red due to their high 
redshifts, there is also a population of red sources where 
the red colours can be attributed to: an intrinsically red 
continuum; significant dust extinction associated with the 
quasar host galaxy; the presence of broad absorption lines; 
or a gravitationally lensed system where the foreground lens- 
ing g alaxy and /or intervening absorption system redden the 
light (iGregg et alfcOOllLacv et al.ll2Q02h . Long wavelengths 
such as the far infra-red and sub-millimetre have for some 
time been used t o track very red, dust enshrouded po pula- 
tions of galaxies JBlain et al.ll2002l : ISmail et al.ll20Q2l ) that 
are otherwise obscured in the optical. Now, the emergence 
of very large area surveys in the near and mid infra-red 
wavelengths such as th e UKIDSS Large Area Survey (LAS; 
iLawrence eraD (l2QQ7l \). VISTA Hemisphere Survey (VHS; 
McMahon et al. in preparation) and Wide Infra-red Survey 
Explorer [WISE] Wright et al. (20i3)), are providing a dif- 
ferent perspective on studies of red populations of galaxies 
and quasars. While these surveys are relatively shallow, the 
huge area of sky surveyed allows detection of rare, luminous 
sources. 

Quasar selection in large surveys has typically relied on 
detecting the ir excess in flux relat ive to stars at short UV 
wavelengths ([Richards et al.ll2002l l. However, quasars also 
displ ay an infra-red excess in the iC-band (Warr en et al.l 
bood ) which allows detection of a complementary sample 
of infra-red bright sources that are rel atively faint at op- 
tical wavelengths ([Maddox et al.l |2008|) . Although numer- 
ous attempts have previously been made at studying the 
infra-red bright quasar population, the prevalence of Galac- 
tic stars with very similar colours to red quasars in near 
infra-red (NIR) surveys has meant that the selection of 
these object s has typically ha d to rely on optical-NIR 
colour cuts (Maddox et alj|2008|) or matching to radio sur- 
veys (JGlikman et al. 2007). A consequence of the heteroge- 
nous selection criteria used to isolate red quasars is that 
the exact fraction of obscured quasars, that are not de- 
tected in optical surveys, remains u nclear. The estimates 
currently ran ge from <20 per cent (IRichards et al. 1 20031 : 
iMaddoxet al. 2008 ) to 60 per cent (iGlikman et ahlboOTf ). 
As dust-obscured quasars are expected to be bright at near 
infra-red wavelengths, a homogenous flux-limited sample se- 
lected in the NIR allows one to estimate their true space 
density without requiring matches to other multi- wavelength 
surveys which inevitably introduce additional selection bi- 
ases on the sample. 

In Hawthorn et al. (2012, In Prep; hereafter Paper I), 
we introduce a technique for isolating stellar Extremely Red 
Objects (EROs) from large area NIR-surveys using sl (J—K) 
colour selection. It is shown that, provided the artifacts that 
dominate colour selected samples of extreme objects are well 
understood, selecting bright EROs with stellar morphologies 
could prove effective in assembling a homogenous popula- 
tion of infra-red selected dusty Type 1 AGN. Spectroscopic 
follow-up of ten sources over ~100 deg^ in Paper I has led to 
confirmation that seven are extremely red Type 1 AGN with 
similar Ha equivalent widths to samples selected in the op- 
tical. The red colours in these quasars imply dust extinction 
values of Ay ^2.3. 



In the current work, we present new spectra from a 
search for stellar EROs over a larger area using the same 
techniques as in Paper I. Taking the samples from both pa- 
pers, we focus on the inferred properties of our infra-red 
bright Type 1 AGN and compare them to quasars selected 
at ot her wavel e ngths ; from the optical to the sub-millimetre. 
Sanders et al.l (| 19881 ) was the first to postulate that sub- 
millimetre luminous galaxies (SMGs) and quasars are dif- 
ferent observational manifestations of the same sources. 
Major merger induced starbursts appear as sub-millimetre 
galaxies and as the dust clears from the decaying star- 
burst, the central nuclear region is revealed as an optically 
bright quasar. Red broad absorption-line (BAL) quasars 
th at are e nshrouded by dust have also been studied e.g. 
bv lEgami e t al. (1996) and hypothesized to correspond to a 
young phase in quasar evolution. A sample of NIR-selected 
red quasars may therefore provide crucial observational evi- 
dence for the link between galaxies and supermassive black 
holes, representing a short-lived phase in the evolution of the 
galaxy where the dust has not yet fully cleared and is still 
obscuring the view of the quasar in the optical wavelengths. 

Throughout this paper we assume a flat concordance 
cosmology with Ho=70 km s~^ Mpc~^, Qm=0.3 and 
^A=0.7. All magnitudes quoted are in the native system 
of the survey being discussed - i.e. AB for SDSS and Vega 
for UKIDSS unless otherwise stated. The AB-Yega conver- 
si ons for both the SP SS and UKIDSS filters can be found 
in lHewett et~all (|2006l ). 



2 DATA 

Red quasar candidates are selected from the UKIDSS LAS, 
one of the largest NIR surveys and the NIR counterpart to 
the Sloan Digital Sky Survey in the northern hemisphere. In 
Paper I, we have presented details of our candidate selection 
algorithm including how to effectively screen the data for the 
many spurious sources. Below, we provide a brief overview 
of some of the key features of our selection but the interested 
reader is referred to Paper I for more information. 



2.1 Photometric Selection 

Red quasar candidates were selected from the UKIDSS- 
LAS Data Release 4 (DR4) by isolating sources with stel- 
lar morphologies and J— Kvega >2.5 [J — Kab >1.5] and 
i^Vega <16.5 [Kab <18.4]. They belong to a class of sources 
known as Extremely Red Objects (EROs). The UKIDSS- 
LAS DR4 footprint is shown in Fig. [T] Two samples were 
defined: 1) 08h<RA<16h and 0°<DEC<15° where we re- 
quire a non-detection in the optical SDSS bands in order for 
the source to be included in our spectroscopic sample, and 2) 
the SDSS southern stripe Stripe82 region at 21 h<RA<00h 
and —1.25°< DEC<1. 25° which benefits from deep co- added 
optical data (jAnnis et al.ll201lh . In the case of the Stripe82 
sample, we use the deep co-add images to select sources that 
are fainter than iy^^=20.5 and redder than {i — K)ab—'2.5 
{{i — iT) Vega =4. 4) as was the case for the North Galactic 
Pole (NGP) part of the Stripe82 sample presented in Paper 
I. We also include four sources that are brighter than this i- 
band limit in order to check whether there are any reddened 
AGN that could potentially be missed by applying the SDSS 
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Figure 1. The UKIDSS-LAS DR4 footprint from which all targets in this paper have been selected. 



i-band selection. Finally, to assess the utility of X-ray de- 
tections for identifying red, luminous A GNs, two objects in - 
cluded in the 2XMM X-ray catalogue (IWatson et al.ll20Q9l ) 
that satisfy Xvega <17.0 [Kab <18.9] were also targeted. 
The objects are the two brightest EROs that are known X- 
ray sources in the survey area. 

The full area spanned by Sample 1 is 770 deg^, cal- 
culated using the publicly available Mangle software 
(ISwanson et al.ll2QQ8l ). There are 43 EROs with i^vega <16.5 
and (J — i^) Vega >2.5. Nine of the 43 sources are undetected 
in the SDSS and therefore are candidate reddened Type 1 
AGN. In addition, ULASJ0908+1042 was targetted despite 
being detected in SDSS as it is also a FIRST radio source. 
Finally, ULASJ1002+0137 was also targetted despite being 
fainter in K as it is an X-ray source. Sample 1 therefore 
comprises 11 red quasar candidates which are summarised 
in Table [H 

The South Galactic Pole (SGP) part of Sample 2 that 
overlaps the SDSS Stripe82 spans an additional 100 deg^. 
There are 14 EROs in this area with ir(Vega)<16.5 and 
(J— i^)vega >2.5. Only three of these are significantly fainter 
than i—20.5 in the deeper SDSS co-add images and redder 
than {i - K)ab=2.5. Two of these, ULASJ2200+0056 and 
ULASJ2224— 0015, were targetted for spectroscopy but the 
third, ULASJ2346— 0038, was unfortunately incorrectly ex- 
cluded from the spectroscopic target list as it had a close 
neighbour in a different frameset and was erroneously as- 
sumed to be a duplicate source. In addition, four other 
bright sources with Kvega <16 that have iAB <=20.5 were 
followed-up in order to assess whether any luminous K- 
band sources that are also detected in SDSS could poten- 
tially be reddened Type 1 AGN. Finally, the fainter source 
ULASJ2343H-0017 was also targetted on account of being 
detected in the X-ray. Therefore Sample 2 comprises eight 
potential red quasar candidates four of which have i<=20.5. 

In Table[T]we provide details of all our red quasar candi- 
dates in both samples including those without spectra. Fig. [2] 
shows the {i — K) distribution of the EROs in Stripe82 with 
deep optical photometry. These include the SGP part of the 



sample presented in Table [T] of this paper as well as the NGP 
part of the sample presented in Paper I. We also plot sev- 
eral well-known red obje cts such as th e dusty ultraluni inous 
galaxies HRIO (IHu fc R idgwav 1994 : IPev etal .11 19991 ) and 



IRAS F10214 (jRowan- Robinson et al 

red lensed quasars such as J0738+2750 and 



19911) 



well-known 
PKS-J0132 



qi^ 

(|Gregg et al.l |2003) and red broad absorption- line qu asars 
e.g. Hawaii-167 (|Cowie et al.lll994l : [Egami et al.lll996l ). As 
can be seen from the figure, our ERO candidates are among 
the reddest infra-red bright objects currently known. As dis- 
cussed in detail in the next section, all the SDSS faint EROs 
are spectroscopically confirmed to be Type 1 AGN whereas 
no redshifts were obtained for the SDSS bright EROs con- 
sistent with identifications as compact red galaxies. 



2.2 Spectroscopic FolloAv-up 

Spectroscopic observations of our red quasar candidates 
were made with the SINFON I integral field spectrograph 
on the VLT (lEisenhauer et al.|[2Q03 : Bonnet et al. 2004 ) be- 
tween 2009 April and July (Period83:383.A-0573(A)) with 
seeing ranging from 0.8-1.4^^ and the sample of 13 red quasar 
candidates that were observed, is presented in Table [1] SIN- 
FONI was used in noAO mode with 0.25^^x0.25^^ pixels cor- 
responding to an 8''x8'' field-of-view with the R=1500 H+K 
grism yielding a resolution of 12 A or 185kms~^ (measured 
from sky lines) and a dispersion of 5 A per pixel. The ob- 
served wavelength range is 1.4-2.45 /xm. The targets were 
offset by ±1.5'' in RA and DEC from the centre of the IFU 
for sky subtraction purposes, however always retaining the 
target in the field of view. The array centre therefore moved 
in a 3x3'' square. The data were reduced using standard 
SINFONI ESOREX pipelines which include extraction, sky 
subtraction, wavelength calibration and flat fielding. 

All spectra were flux- calibrated using standard telluric 
stars observed at similar airmass. The target spectrum is di- 
vided by the spectrum of the calibration star and multiplied 
by the appropriate blackbody, depending on the spectral 
type of the star, in order to get a relative flux- calibrated 
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Table 1. Sample of red quasar candidates from the UKIDSS-LAS DR4. Two samples are constructed at 08h<RA<161i and overlapping 
the SGP part of the SDSS southern stripe Stripe82. Slightly different selection criteria were used to select sources in these two samples 
as detailed in Section 12.11 Both sources that were spectroscopically followed-up using SINFONI as well as sources with no spectra are 
included for completeness. The NGP part of the Stripe82 sample was presented in Paper I as well as two of the quasars in Sample 1 that 
overlap the UKIDSS-LAS Early Data Release area. 



Name 



RA 



DEC 



Kv. 



ega 



lAS 



Redshift 



Notes 



Sample 1 08H<RA<16H - 11 objects 



ULASJ0908+1042 09:08:06.81 +10:42:39.5 15.13 

ULASJ0946+0016 09:46:28.66 +00:16:25.8 16.28 

ULASJ1002+0137^ 10:02:11.29 +01:37:06.8 17.00 

ULASJ1004+0026" 10:04:43.68 +00:26:09.8 15.51 

ULASJ1200+0423 12:00:59.45 +04:23:38.4 16.30 

ULASJ1234+0907 12:34:27.52 +09:07:54.2 16.15 

ULASJ1319+0009 13:19:10.66 +00:09:56.1 16.04 

ULASJ1320+0948 13:20:47.95 +09:08:05.8 16.46 

ULASJ1422+1023 14:22:39.13 +10:23:34.2 16.46 

ULASJ1455+1230 14:55:21.00 +12:30:08.6 16.34 

ULASJ1539+0557 15:39:10.16 +05:57:49.7 15.85 



21.00 - FIRST radio source 

- no spectrum 

1.595 Faint K-band, 2XMM X-ray source 



no spectrum 
2.503 

no spectrum 

1.460 
2.658 



UKIDSS-EDR Paper I 
FIRST radio source 
UKIDSS-EDR Paper I 



Sample 2 Stripe82 SGP 21H<RA<24H - 8 objects 



ULASJ2148-0011 21:48:25.35 -00:11:32.1 15.49 

ULASJ2150-0055 21:50:48.27 -00:55:45.8 15.13 

ULASJ2200+0056 22:00:24.87 +00:56:04.8 15.22 

ULASJ2219+0036 22:19:30.44 +00:36:26.4 15.61 

ULASJ2224-0015 22:24:09.41 -00:15:24.1 16.07 

ULASJ2251+0047 22:51:21.84 +00:47:33.4 15.65 

ULASJ2343+0017 23:43:39.68 +00:17:56.0 16.68 

ULASJ2346-0038 23:46:58.41 -00:38:06.4 15.31 



20.41 - Bright SDSS 

20.53 - On SDSS selection hmit 

21.00 2.541 Faint SDSS 

19.06 - Bright SDSS 

21.77 2.223 Faint SDSS 

20.08 - Bright SDSS 

22.00 - Faint K-band, 2XMM X-ray source 

21.60 no spectrum Faint SDSS 
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Figure 2. i — K vs iab colour-magnitude diagram showing the 
EROs from our survey overlapping the deep S DSS Stripe82 co- 
add d ata. We also show the SDSS DR7 quasars (|Schneider et al.l 
l2010h as the grey p oints, Luminous Re d Galaxies (LRGs) from 
the 2SLAQ survey (|Cannon et al.l 120061 ) as the small red points 
as well as various well-known red sources in green. Our selection 
criteria that divides SDSS-bright and faint EROs is shown by the 
dashed lines. The SDSS-bright EROs with i<20.5 (blue) have no 
lines in the NIR spectra consistent with these being LRGs whereas 
the SDSS-faint EROs (red) are all spectroscopically confirmed to 
be reddened Type 1 AGN. 
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Figure 4. SINFONI spectrum of our brightest unconfirmed tar- 
get UL AS J 1004+0026 which is unlikely to be an elliptical or low- 
redshift radio galaxy. The atmospheric transmission spectrum 
from VLT-ISAAC is also shown for reference at the top of the 
panel and demonstrates the wavelength region (1.80-1.97/im cor- 
responding to Ha at 1.75< z <2) where no lines can be observed 
due to atmospheric opacity. 



spectrum. Due to the absence of absolute spectrophotomet- 
ric standards in the NIR, the relative flux- calibrated spectra 
are then normalised to the UKIDSS K-hdJid magnitudes of 
the targets in order to provide an absolute flux calibration. 
Three of the 13 quasar candidates observed show broad 
Ha emission in the iT-band and two have Ha emission in the 
if -band. We fit multiple Gaussian components to the Ha 
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Figure 3. Near infra-red spectra of all 12 confirmed reddened quasars with redshifts . The top two rows show VLT-SINFONI spectra 
apart from ULASJ1539+0557 which was observed using NIRSPEC on Keck. The bottom two rows show Gemini-GNIRS spectra. In the 
case of the GNIRS spectra, the H and X-band orders are plotted together in these figures but problems with extracting the spectra at 
the edges of each order mean that the relative fiux calibration between the two orders is unlikely to be accurate. All spectra have been 
smoothed by a 3-pixel boxcar (0.0015 /im) for purposes of presentation. The red dot-dashed lines show the model fits derived for each 
quasar by fitting to the broadband colours. The model SEDs are derived by fitting to the broadband colours of the quasars with a simple 
scaling of Ho; equivalent width to match the observations. The Ha line profiles in the models and observations are therefore not expected 
to match. The region of low atmospheric transmission between the H- and X-bands is shown as the grey rectangular area. 



line profile and redshifts are determined from the centroid 
of the intermediate component where present and from the 
centroid of the single Gaussian where only one component is 
required to fit the data. Redshifts are summarised in Table 
(TJ In Fig. [3] we show spectra of the five new red quasars 
with redshifts as well as the seven confirmed quasars already 
presented in Paper I. 



Out of the eight sources for which no redshifts were 
obtained, two have associated FIRST counterparts so these 
could be lower redshift radio galaxies with compact mor- 
phologies. One of the X-ray sources is at the faint-end of 
our sample and below the nominal flux limit. However, the 
other X-ray source, ULASJ1002+0137, which also lies be- 
low the nominal flux limit of the survey was spectroscopi- 
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Table 2. Photometric redshifts for the four SDSS Stripe82 EROs 
with i <=20.5 and for which no hnes were identified in the NIR 
spectra. 



Name Zp 

ULASJ2148-0011 0.97±^:^^ 

ULASJ2150-0055 1.30±0.08 

ULASJ2219+0036 2.71±0.06 



ULASJ2251+0047 0.85±0.05 



cally confirmed to be a broad-line quasar at z = 1.595. This 
source is in the COSMOS field and therefore has ancillary 
multi- wavelength data. We discuss it in detail in § 13.31 

In Sample 2, all the sources that had iab >20.5, 
were spectroscopically confirmed to be reddened Type 1 
AGN. No redshifts were obtained for the four sources 
with iAB <=20.5, consistent with these being com- 
pact elliptical galaxies. We note that three out of the 
four sources - ULASJ2148-0011, ULASJ2150-0055 and 
ULASJ2251+0047 - are indeed classified as galaxies in the 
SDSS imaging survey. We u se the publicly availa ble photo- 
metric redshift code Eazy (|Brammer et al.ll2008l l to calcu- 
late photometric redshifts for these four SDSS-bright EROs 
using an R-hand luminosity function prior. The photometric 
redshifts of the four sources are summarised in Table [2] 

The sources without redshifts which are also undetected 
in the optical filters could have several interpretations. The 
spectrum for the brightest of these, ULASJ1004+0026, is 
shown in Fig. |4] along with the sky transmission spectrum. 
If such unconfirmed candidates are quasars, as suggested 
by their compact morphologies and NIR colours, they may 
be at redshifts between 1.75-2.00 which corresponds to Ha 
lying in the region of strong atmospheric transmission be- 
tween the H and K bands. Alternatively, there is a class 
of ra re broad absorp t ion line quasars known as FeLOB- 
ALs (IVoit et al.1 Il993l : iBecker et al.] Il997l . l2000l : iHall et al.1 
I2OO2I I where strong absorption due to iron at relatively low 
ionisation potentials means there is virtually no flux below 
A ^ 2800 A. It has been hypothesised that FeLOBALs may 
represent a transition population between ultra-luminous 
infra-red galaxies and quasars with some of these sources 
having prodigiously high star form ation rates, as inf erred 
from their mid infra-red properties ([Farrah et al.ll2007l V Al- 
though they are very rare, they would be preferentially se- 
lected in N IR-selected sa r aples due to their very red colours 
- e.g. as in lUrrutia et al.l (|2009|). If some of our unidentified 
objects do indeed correspond to this class of objects, this 
would place them at z ^3-4 in order to reproduce the very 
red colours. This is consistent wit h the fraction of BALs 
increasing at these higher redshifts (| Allen et al.ll201lh . Op- 
tical spectroscopy would help confirm the identity of these 



3 RESULTS 

The spectra presented above combined with the sample in 
Paper I total 12 confirmed Type 1 AGN with very red 
colours and redshifts between 1.4-2.7. We now turn to a 
detailed assessment of their properties. 



3.1 Model SEDs 

We begin by fitting model Type 1 AGN SEDs to the 
available SDSS+UKIDSS broadband photometry for all our 
quasars. The photometric data is chosen for the fitting as 
it covers a larger wavelength range than our near infra- 
red spectra and includes flux points at bluer wavelengths 
which are more sensitive to dust extinction. The unred- 
dened 'standard' qua sar model is described in Section 2.4 of 
iMaddoxetaD (|2008l ). The model SED reproduces the SDSS 
ugriz and UKIDSS YJHK broadband colours of the ma- 
jority of quasars, selected in the SDSS, to high accuracy 
over the redshift range 0.2 ^ z ^ 5. The aim of fitting 
these model SEDs to the quasar photometry, is to infer the 
amount of dust extinction in each quasar from the observed 
colours. Although other factors such as lensing and variabil- 
ity may also contribute to the red colours, these are unlikely 
to be significant in our sample. All our quasars lie at z>1.5 
and are very luminous so the host galaxy contribution to the 
broadband light is likely to be insignificant. While variability 
can certainly affect the spectral shape of quasars, the effect 
is more pronounced in the rest-frame UV bands and lower 
luminosity objects (IHook et al.lll994ll and mo re often seen 
in radio- loud populations (|Helfand et al.ll200lh . Only two of 
our 12 quasars, ULASJ0016-0038 and ULASJ0141+0101 
are radio- loud. We discuss the possibility of lensing in detail 
in § 14.11 where we conclude that it is unlikely to affect the 
majority of our sample. We therefore proceed by assuming 
that the red NIR colours seen in our quasars is due to dust 
extinction. 

Reddened model SEDs were therefore created for each 
quasar, by including the effect of an extinction curve with a 
specified E{B — V). T he form of the extincti on curve used is 
very similar to that of lGallerani et al.l (|2010|) but in the rest- 
frame wavelength interval of interest here, 2500-10 000 A, 
the differences between commonly used extinction curves 
(e.g. LMC, SMC and Milky Way) are small and the results 
are not dependent on the choice of extinction curve. The 
facility to adjust the emission line equivalent widths is also 
incorporated in the models and object-specific model SEDs 
were generated by matching the measured Ha equivalent 
widths. We stress however, that we do not attempt to match 
the Ha line profiles in the model and observed spectra apart 
from this simple scaling of the equivalent widths. In this way, 
we can disentangle the effects of Ha equivalent width and 
dust extinction in the observed (J — K) colours. 

In a few cases, the {J — K) and {H — K) colours are 
marginally inconsistent with a single reddening value so we 
quote the best-fit extinction values needed to reproduce each 
of the observed colours. However, the inferred E{B — V) 
values are always similar and the conclusions of the paper 
do not depend critically on which values are adopted. For 
the rest of the paper, we use the extinction values derived 
by matching to the (J — K) colours so as to provide a larger 
baseline in terms of wavelength for the extinction estimates. 
An Rv=3.1, consistent with the dust in our own Milky Way, 
is assumed throughoutij. 

In Table [3] we quote the dust extinction. Ay inferred 
from the model fits for our combined sample of reddened 



^ Considering possible variations in the dust-law in distant 
quasars is beyond the scope of this work. 



© 0000 RAS, MNRAS 000, 000-000 



7 



quasars from Paper I and the current work. We also give 
the factor by which the Ha equivalent widths in the model 
are scaled in order to rep r oduce the observations. Previous 
studies by IClikman et al.l (|2007|) did not account for the ef- 
fect of the Ha equivalent widths on the {J — K) colours. 
Fig. [3] shows these best-fit model SEDs overlaid on the ob- 
served rest-frame optical spectra for all our spectroscopically 
confirmed red quasars. Note that the model fits were derived 
from the broadband photometry only with a simple scaling 
adopted to match the observed Ha equivalent widths so the 
match between these model SEDs and the observed contin- 
uum flux in the spectra is, in most cases, very good. In some 
cases - e.g. ULASJ0221— 0019 - the models do not match the 
continuum in the observed spectra at the edges of the spec- 
tral orders due to problems with extracting the spectra at 
the edges. We see from Table [3] however that the match be- 
tween the observed and model colours is always very good. 
The median difference in the model and observed (H — K) 
colours is 0.005. 

From Table [3] it can immediately be seen that most of 
the red quasars have similar Ha equivalent widths to the 
standard model spectrum constructed from the optically se- 
lected quasar population, with the strongest Ha line having 
an equivalent width of only a factor of two larger. However, 
we find that significant dust extinction is required to repro- 
duce the continuum colours, with a mean dust extinction of 
Av --2.5 for ah our sources. ULASJ1234+0907, our reddest 
object, has a dust extinction of Ay=6.0 inferred from its 
[H — K) colour and is one of the reddest quasars currently 
known. 



3.2 Line Widths, Bolometric Luminosities &; 
Black Hole Masses 

In this section we assess the emission line properties of our 12 
red quasars. In order to do so, we fit multiple Gaussian com- 
ponents to the Ha line profiles of all our quasars with spec- 
troscopic redshifts after blueshifting to the rest-frame. Note, 
these Gaussian fits to the observed Ha lines are distinct from 
the model SED fits done to the broadband photometry in § 
l3.1l to infer a dust extinction. The effect of this dust extinc- 
tion on the full width half maximum (FWHM) of the Ha 
line is expected to be very small. The Gaussian fits can be 
seen in Fig. [5] and Fig. [6] for the new sample presented in 
this paper and the Paper I sample respectively. The derived 
line- widths are summarised in Table |4l Errors are derived 
from a least-squares fit of the Gaussian profiles to the data. 
In most cases, both intermediate and broad components are 
needed to fit the line profiles adequately. The presence of 
a very broad Ha component with FWHM >5000kms~^ in 
these quasars, confirms that we are able to at least partially 
see through into the broad line region located close to the 
black-hole accretion disk. The broad components in our sam- 
ple typically have FWHM=6000kms~^ while the intermedi- 
ate component has an average FWHM=2230kms~^. Our 
red quasars are therefore likely to be canonical Type 1 AGN 
rath er than Type 2 AGN that are ob scured by a dusty torus 
(e.g. iMartinez-Sansigre et al.l ((20051)). In the latter case we 
would not be able to detect the broad component for most 
lines of sight. We conclude that the large dust extinction val- 
ues must be arising in the host galaxy of the quasar rather 
than from the molecular torus. 



The Ha line properties can be used to infer black- 
hole masses, bolometric luminosities as well as Eddington 
fractions usin^ the scaling re lations between Ha and H/3 
linewidth (IGreene k, Hdl2005l ) from reverberation mapping 
([Kaspi et al.ll2007l V We correct the line for the instrumental 
broadening of 185kms~^ by subtracting this in quadrature 
from the observed linewidth and use the equations below: 



FWHMh/3 = (1.07±0.07) X 10^ ( ^" ) kms"' 

V lO^kms / 

(1) 

The black hole mass is then determined from the H/3 FWHM 
using: 



Mbh/Mct 



10^ 



/ FWHMh/3 
VlOOOkms-^ 



lO^^ergs- 



(2) 



where L5100 is estimated using the best-fit model SEDs of 
each quasar. We assum e that the bolometric luminos ity is 
given by Lboi^TLsioo ([Vestergaard k, PetersonI |2006|) and 
use these to calculate the Eddington fraction. The mass ac- 
cretion rate is given by: 



^=0.18^ 
dt e 



1046 



ergs 



Moyr- 



(3) 



where the accretion efficiency, e is assumed to be 10 per cent 
and the factor of 0.18 results from the units being used in 
EqO The derived luminosities, masses, Eddington fractions 
and accretion rates are summarised in Table [5] The quoted 
errors are derived by propagating the fit errors to the Ha 
linewidths. 

In addition to these fit-errors, there are also systematic 
errors affecting the black-hole mass estimates, that should 
be discussed. We therefore highlight several caveats relevant 
to the black- hole mass estimates in Tabled which may re- 
sult in large systematic uncertainties. In all cases, black-hole 
masses were estimated using the FWHM of a single Gaussian 
fit to the Ha line profile. However, using the broad compo- 
nent only to estimate the black-hole masses would result in 
an increase in the black-hole masses and correspondingly a 
decrease in the Eddington fractions by up to a factor of ^2. 
Furthermore, we note that several of our Gaussian line pro- 
files in Fig.[5]and Fig.[6]are highly asymmetric showing large 
velocity shifts between the different components. 

3.2.1 Evidence for Outflows: 

In most cases, the broad component shows a clear blueshift 
relative to the intermediate component. This blueshift is 
most pronounced in ULASJ0016-0038 and ULASJ0041-0021 
where the broad component has a velocity of -1400 kms~^ 
and -800kms~^ relative to the intermediate component in 
the two respective sources. In ULASJ1455+1230 however, 
the broad component is redshifted by +760kms~^ rela- 
tive to the intermediate component. Outflowing material 
in a bipolar flow can produce both blue and red-shifted 
emission so these velocity shifts indicate the presence of 
strong outflows that are affecting the broad line region 
of the quasar. Signatures of outflows, while not generally 
seen in the UV-luminous quasar population, are not uncom- 
mon in samples of red BAL quasars where they have been 
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Table 3. Summary of Model Fits for Spectroscopically Confirmed Red Quasars. 



Name 


Redshift 


{J - KUs 


{J — ^)model 


{H - X)obs 


{H — X)model 


Ay 


Ha EW Scaling 


ULASJ0016-0038 


2.194 


1.607 


1.702 


0.871 


0.882 


1.7 


2.1 


ULASJ0041-0021 


2.517 


2.171 


1.821 


0.851 


0.869 


2.1 


0.7 




2.517 


2.171 


2.167 


0.851 


1.023 


2.5 


0.7 


ULASJ0141+0101 


2.562 


1.767 


1.752 


0.893 


0.893 


1.8 


1.2 


ULASJ0144-0114 


2.505 


1.852 


1.852 


1.009 


0.897 


2.1 


1.0 




2.505 


1.852 


2.099 


1.009 


1.007 


2.4 


1.0 


ULASJ0144+0036 


2.283 


1.902 


1.496 


0.726 


0.734 


1.7 


0.9 




2.283 


1.902 


1.901 


0.726 


0.913 


2.3 


0.9 


ULASJ0221-0019 


2.247 


- 


1.690 


0.854 


0.856 


1.8 


1.6 


ULASJ1539+0557 


2.658 


- 


3.290 


1.440 


1.450 


4.0 


0.3 


ULASJ1002+0137 


1.595 


- 


1.657 


0.551 


0.550 


3.2 


0.8 


ULASJ1234+0907 


2.503 


- 


5.065 


2.321 


2.331 


6.0 


1.3 


ULASJ1455+1230 


1.460 


2.212 


1.729 


0.570 


0.571 


3.4 


0.9 




1.460 


2.212 


2.219 


0.570 


0.790 


4.3 


0.9 


ULASJ2200+0056 


2.541 


1.693 


1.617 


0.798 


0.815 


1.7 


1.3 


ULASJ2224-0015 


2.223 


1.747 


1.798 


0.910 


0.912 


1.9 


1.7 



Table 4. Ha FWHM derived by fitting both a single Gaussian as well as two (broad and an intermediate components) Gaussian to the 
Ha line profiles. A narrow Ha component was not required in any of the sources. 



Name 



ULASJ0016-0038 
ULASJ0041-0021 
ULASJ0141+0101 
ULASJ0144+0036 
ULASJ0144-0114 
ULASJ0221-0019 
ULASJ1539+0557 
ULASJ1002+0137 
ULASJ1234+0907 
ULASJ1455+1230 
ULASJ2200+0056 
ULASJ2224-0015 



FWHM^^^^g^'^/kms-^ FWHM^"°^^/kms-i FWHM^^^^/kms" 



5400±400 
5500±500 
5900±300 
6000±200 
4700±400 
3700±200 
4700±400 
3200±300 
7300±600 
2900±200 
3300±200 
3100±100 



7900±200 
8700±600 
7300±300 
6400±200 
5400±500 
6000±100 
5900±400 

8100±600 
6400±300 
6200±100 
6700±200 



2620±70 
2800±200 
3200±100 

1410±40 
2800±200 

1950±30 
2200±200 

2600±200 
2000±100 
1960±40 
2350±60 



sho wn to quench star formation within their host galax- 
ies (|Farrah et al.1 l2012l l . These flows can broaden the Ha 
linewidths leading to an over- estimation of the virial black- 
hole masses and consequently an under-estimate of the Ed- 
dington fractions. We do not correct for this effect in our 
black-hole mass estimates but caution that excluding the 
blue/red asymmetric component from the Ha line fits in 
quasars showing evidence for outflows, results in a decrease 
in the estimated black-hole masses by a factor of ^2-3. 



to the kn own lensed reddened quasar PKS01 32 at similar 
redshifts JHall et alJ l2002l : iGregg et al1l2002l ). Our ongo- 
ing follow-up observations of ULASJ 1234+0907 at long 
wavelengths will enable a detailed characterisation of its 
host galaxy properties, shedding light on the reason for 
the extreme luminosity and dust extinction. Although we 
cannot exclude lensing as a possible reason for the extreme 
luminosity of this source we discuss why this is statistically 
improbable in Section (ITTl 



3.3 Notes on Individual Sources 

Two of the sources in our reddened AGN sample are 
particularly interesting and merit further discussion: 

ULASJ1234+0907 : ULASJ1234+0907 is our reddest 
quasar with {H — K)=2.3 and an inferred dust extinction of 
Aa/=6.0 derived from this colour. The extinction corrected 
bolometric luminosity for this source is Lboi >10'^^erg s~^ 
and the inferred black-hole mass is >10^^Mq from the very 
broad Ra line seen in Fig. [3 This makes ULASJ1234+0907 
one of the most bolometrically luminous and massive 
black-holes known at z ~2. As discussed in detail in 
Section [J the extinction-corrected absolute magnitude for 
ULASJ1234+0907 is M^ - -31 — 32 which is comparable 



ULASJ1002+0137 : ULASJ1002+0137 is in the COS- 
MOS field and an d has previously been studied by 
iMainieri et al.1 (|2007l ) in their analysis of X-ray AGN from 
the XMM-Newton Wide Field Survey. The o ptical coun- 
terpar t to the X-ray source is presented in iBrusa et al.l 
(|201Q| ). Although a redshift of z =0.784 was initially de- 
rived bv IMainieri et al.1 (|2007l ) and lTrump et all (l2009|). the 
optical spectroscopic redshift of z =1.592 in iBrusa et al.l 
(2010) agrees well with our Ha re dshift and the photomet- 
ric redshift of I Sal vat o et al.l ((20091). The X-ray spectrum is 
best fit by an absorbed power-law and an Fe Ka line at 
z =0.784. The multi-wavelength analysis and optical spec- 
troscopy prod uces a classification as a Narrow-Line/Type 
2 AGN (Brusaet al.1 [2010- ). Our measured Ra FWHM of 
^3000 kms~^ would suggest that the broad- line region is at 
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Figure 5. Gaussian fits to Ha line for the five new quasars. All observed spectra have been smoothed by a 5-pixel boxcar for purposes of 
presentation. A broad as well as an intermediate component are needed to fit the line profiles adequately for most of the sources. These 
individual components are shown as the dashed lines whereas the solid lines show the sum of the two components. We note that in the 
case of ULASJ1002+0137 a single-component fit, with FWHM=3200±300kms~\ produced an equally good x^ • The residuals from the 
fits are also shown at the bottom of each panel. 



least partially visible. However, this FWHM is significantly 
smaller than the broad components for the rest of the sam- 
ple and ULASJ1002+0137 is the only source where a single 
Gaussian provides an adequate fit. If the broad- line region 
is partially obscured by a dusty torus, the large dust ex- 
tinction of Av=3.2 inferred for this source would mean that 
the broad-line flux is expected to be more suppressed in the 



rest-frame ultraviolet than in the optical. This may have 
prevented this source from being classified as a broad-line 
AGN from the optical spectrum. 

We should also consider the possibility that this is a 
lensed quasar with the source at z =1.595 and the lens at 
z =0.784. In Fig.Elwe show a 15x15'' COSMOS BiK com- 
posite around the source at a resolution of 0.15'' per pixel 
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Table 5. Bolometric luminosities, virial black-hole masses, Eddington fractions and mass accretion rates derived for the 12 reddened 
quasars. 



Name 


logio(Lboi/erg s-^) 


logio(MBH/M0) 


Lboi/LEdd 


dm/dt (Mq yr 


-') 


ULASJ0016-0038 


46.9 


9.32±0.06 


0.20±0.05 


13 




ULASJ0041-0021 


47.3 


9.56±0.08 


0.30±0.01 


37 




ULASJ0141+0101 


46.9 


9.41±0.04 


0.17±0.02 


14 




ULASJ0144-0114 


47.3 


9.47±0.07 


0.37±0.02 


36 




ULASJ0144+0036 


47.1 


9.53±0.03 


0.19±0.01 


22 




ULASJ0221-0019 


46.8 


8.96±0.05 


0.40±0.02 


12 




ULASJ1539+0557 


48.2 


9.87±0.07 


1.07±0.06 


264 




ULASJ1002+0137 


46.8 


8.83±0.08 


0.50±0.03 


11 




ULASJ1234+0907 


48.5 


10.43±0.07 


0.65±0.04 


580 




ULASJ1455+1230 


47.2 


8.94±0.06 


1.04±0.05 


30 




ULASJ2200+0056 


47.4 


9.18±0.05 


0.99±0.04 


50 




ULASJ2224-0015 


47.0 


8.93±0.03 


0.68±0.02 


19 






Figure 7. COSMOS B, i, K composite 15x15^^ cut-out around 
ULASJ1002+0137. The quasar is marked with the circle. An ex- 
cess of blue galaxies is seen around the source and the three 
sources located southwest of the quasar could be a lensed arc. 



(jMcCracken et al.ll20ld V An excess of blue sources is seen 
around the red quasar and the three blue sources to the 
southwest of ULASJ1002+0137 could be a lensed arc al- 
though the suggestion is highly tentative. We discuss the 
possibility of lensing in our sample further in Section 14.11 



4 DISCUSSION 

4.1 Lensing in the Reddened Quasar Population 

It has previously been demonstrated that lensing may be 
responsible for reddening the opt ical and infra-red colours 
of quasars ([Malhotra et al.l 119971 ). There are various ex- 
planations for why this may be the case. Firstly, dust 
from a lens galaxy could extinct the light coming from 
the quasar as in the cas e of the lensed quasar PKS1830- 
211 (jCourbin et al.l [l998l V Alternatively, there could be a 
contribution to the quasar colours from galactic starlight 
within the lensing galaxy. There are various examples of 
lensed quasars in the literature which were selected based on 
their infra-red excess. These include JO 134-0931 oth erwise 
known as PKS Q132 (IWinn et al. 1 l2QQ2l: iGregg et al. 



20021 : 



iHall et aLlliooi ). F IRSTJ1 004+1229 JLacv et al.ll2002l )md 
SDSSJ1313+5151 (lOfeket al. 20071 ). In Fig.[8]we compare 
the colours and magnitudes of our reddened quasars to some 
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Figure 8. (J — K) versus K (Vega) colour-magnitude diagram 
comparing our sample of reddened quasars to some well-known 
lensed reddened quasars at similar redshifts as well as lower red- 
shift LRGs, representative of the likely lens. Our sample of red- 
dened quasars are always fainter and almost always redder than 
the known lensed reddened quasars and significantly redder than 
the LRG lenses. 



of these known lensed quasars, selected based on their infra- 
red excess. In the same plot we also show the distribution of 
LRGs from the 2SLAQ survey at 0.4< z <0.8 (ICannon et al.1 
[2006) which have been matched to the UKIDSS-LAS using 
a matching radius of 1.5" . These LRGs may be considered 
to be representative of the likely lens population. 

In the previous section we presented some very tentative 
evidence for ULASJ1002+0137 being a lensed object based 
on the two different redshifts measured and the excess of 
blue galaxies in the deep COSMOS images. It is interesting 
to consider the likelihood of lensing in our population of 
NIR selected red quasars. Based on Fig. [8] it can be seen 
that our reddened quasars are both apparently fainter and 
redder in colour than most of the known lensed red quasars. 
PKS 1830-211 has a very red observed colour on account of 
lying in the Galactic plane but this is not generally true 
of the other lensed objects. Lensing galaxies are typically 
LRGs at z <1 in which case they should be detectable in 
the SDSS i and z bands as well as the bluer UKIDSS bands. 
The reddened quasars that are undetected in the SDSS as 
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well as the UKIDSS Y and J-bands include our reddest and 
intrinsically brightest source ULASJ1234+0907. Even if the 
lensing galaxy had a redshift as high as z :^1, we would 
expect to be able to detect it in the UKIDSS Y and J-bands. 
We conclude that although we cannot rule out the pos- 
sibility of some of the reddened quasars in our sample being 
lensed, in general this is unlikely to be the case for the bulk 
of the population and dust extinction within the quasar host 
is a more plausible explanation for the extreme colours ob- 
served in these sources. 



4.2 Reddened Quasars: A Transition Population 
Between the Starburst &; Optical AGN 
Phases? 

Our study of extremely red near infra-red selected bright 
quasars has unearthed a population of broad-line AGN at 
z r^2 corresponding to the main epoch of galaxy forma- 



tion with significant dust extinction of Av 



> 



2mags and 



comparable to that measured in the sub-millimetre galaxy 
population ([Takata et al.ll2006l V We have already ruled out 
lensing as a possible explanation for the extreme colours in 
the bulk of this population. The very large observed widths 
(^5000 km/s) of the Balmer emission lines support the hy- 
pothesis that the reddening is produced by dust within the 
host galaxy rather than a molecular torus which would ob- 
scure our view of the AGN broad-line region. An explana- 
tion that does not rely on viewing orientation is supported 
by the fact that the bolometric luminosities are typically 
very high and comparable to some of the most luminous 
optically selected QSOs. The inferred mass accretion rates 
are therefore also high. The fuel for this black-hole accre- 
tion likely comes from intense star formation within the host 
galaxy consistent with the high dust extinctions being pro- 
duced by a major starburst event. Highly star- forming, dust- 
obscured, ultra- luminous infra-red galaxies (ULIRGs) and 
quasars are thought to be different observational manifesta- 
tions of the same phenomena (Sanders et al. 1988 ). Major- 
merge r induced starbursts appear as ULI RGs or SMGs at 
z -2 (JBlain et al.ll2002l : ISmail et al.]l2002l ) and as the dust 
clears from the decaying starburst, the central nuclear re- 
gion is revealed as an optically bright quasar. Our sample of 
infra-red bright reddened quasars may therefore represent a 
transition phase between the starburst and optical quasar 
phases where the dust has not yet fully cleared prevent- 
ing the sources from being selected as ultraviolet-luminous 
quasars. This blowout phase is also evidenced by the pres- 
ence of strong outflows that seem to be affecting the Ha line 
profiles in some of our quasars. 

We now compare our sample properties to those of 
AGN samples selected at other wavelengths. These in- 
clude the optically sel ected quasar population from SDSS 
DR7 (IShen et al.l 201lh , s ub-millimetre bright q uasars from 
ICoppin et al.1 (l2Q08l ) and lOrellana etHI (|201lh and a sub- 
set of sub- millimetre galaxie s hosting AGN tha t were se- 
lected from X-ray stu dies (lAlexander et al.ll2Q08l V The sub- 
millimetre quasars in lOrellana et al.l (1201 ih are all optically 
bright blue quasars whereas those in ICoppin et al.l (|20Q8l l 
include both standard blue optical quasars as well as redder 
X-ray absorbed quasars and quasars detected in blank-field 
sub-millimetre surveys. In Fig. [9] we plot the black-hole 
masses as a function of redshift as well as bolometric lumi- 



nosity for all these samples. The AGN bolometric luminosity 
for the SMGs was est imated by scaling the X-ray luminosi- 
ties by a factor of 35 ([Alexander et al.ll2008l V Note, that the 
scaling relations we use implicitly assume that the black- 
hole mass correlates with the bolometric luminosity so the 
trend seen in Fig. [9] should be interpreted with caution. 

Fig.[9]shows that our sample of reddened quasars consti- 
tute some of the most massive AGN at z ^2 w ith black-hole 
masses com parable to the quasars studied bv lCoppin et al.l 
(200g) and lOrellana et al.l ([2011) and an order of magni- 
tude larger than SMGs hosting AGN at similar redshifts. 
We n ote in particu l ar tha t the most massive black-holes 
in the ICoppin et al.l (|2008|) sample are the optically bright 
subset of their population and the X-ray absorbed and sub- 
millimetre detected sources in their sample have relatively 
modest black-hole masses. Our quasars are therefore as mas- 
sive as the most luminous optically selected quasars but still 
have significant amounts of dust. These results are consistent 
with a scenario where the growth of the black-hole follows 
the starburst phase. In other words, the quasars found in 
blank sub-millimetre surveys and the SMGs hosting AGN 
have not yet had time to significantly grow their black holes 
and represent an earlier evolutionary stage in the lifecycle 
of these sources. 

Recent Iv lHickox et al.l (|2012i ) used clustering arguments 
to provide indirect evidence that powerful starbursts and op- 
tical quasars occur in the same systems. Our ongoing follow- 
up of the reddened quasars at submillimetre wavelengths 
would provide direct evidence for this link between star- 
bursts and AGN. If the reddened quasars are indeed transi- 
tioning from starbursts to UV-bright quasars, the expecta- 
tion is that the host galaxies of our population are likely to 
be more luminous at far infra-red and submillimetre wave- 
lengths than the hosts of the optical QSO population (e.g. 
[isaak et all (l2002l ) and IPriddev et al.1 (|2003l \). If we assume 
that most of the mass in quasars is assembled at 2; >= 2 
(jMcLure fc Dunlop||2004l l consistent with downsizing in the 
quasar population, then the reddened quasars are likely be- 
ing observed at a younger stage than UV luminous quasars of 
similar mass, which have presumably already been through 
an obscured phase associated with the bulk of their mass 
assembly. 

New large-area far infra-red to millimetre surveys such 
as Herschel- ATLAS and surveys with the South Pole Tele- 
scope, could offer the opportunity to directly detect the host 
galaxies of this rare population due to surveying much larger 
areas than covered by previous surveys of the SMG popula- 
tion. 



4.3 Obscured Fraction of Type 1 AGN 

We now attempt to derive some constraints on the space 
density of obscured Type 1 AGN in order to address whether 
these optically obscured sources constitute a significant frac- 
tion of the total Type 1 AGN population. Previous studies of 
highly obscured broa d- line quasars have re li ed on mat chin, 
to rad io surveys (e.g. I Webster et al.l (| 19951 ): iGlikman et al 
(|2007|)) so it is interesting to consider whether similar con- 
clusions about the obscured fraction are reached by study- 
ing the more numerous radio-quiet quasars. Obscured, in the 
context of the following discussion, is taken to mean quasars 
with Av > 2. 
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Figur e 9. Left: Redshift versus black-hole mass of our reddened quasars compared to sub-millimet re bright quasars from I Coppin et aU 
(|2008l ) and lOrellana et al] (|201lh as well as X-ray selected sub-millimetre galaxies hosting AGN from [Alexander et alj (|20Q8l V Our sample 
of reddened quasars constitute some of the most massive black holes at z ^2. Right: Bol ometric luminosit y versus black-hole mass for our 
sample of reddened q uasars compared to the optically selected qua sars from the SPSS (|Shen et aLll201ll ), sub-millimetre bright quasars 
(|Orellana et al.lfeoilh and sub-millimetre galaxies hosting AGN ([Alexander et al.l 120081 ). The dashed line corresponds to Eddington 
accretion. Note, the scaling relations used to derive the black-hole masses and bolometric luminosities means the two quantities are 
necessarily correlated. 



In Section [3. II we fitted model SEDs to the broadband 
colours of our reddened quasars in order to infer a dust ex- 
tinction for each quasar. Using this dust extinction value 
and our reddening law (Section 13. ip , we correct each model 
SED in order to derive an unreddened SED for each quasar. 
The left-panel of Fig. [10] shows an example of this proce- 
dure for our reddest quasar, ULASJ1234+0907. It can im- 
mediately be seen that the extinction in the SDSS optical 
bands is significant; an Av=6.0 corresponds to an extinction 
of 17.5 magnitudes at rest-frame 2250 A which makes the 
quasar undetectable in even deep ground-based optical sur- 
veys. Correcting for the extinction makes ULASJ1234+0907 
one of the brightest quasars known at these redshifts with 
an absolute magnitude of Mi:=-31.8. In Fig. [10] we also plot a 
model with Av=4..96 for ULASJ1234+0907 which produces 
a slightly better match to the spectrum continuum. Adopt- 
ing this lower extinction value rather than the one inferred 
from fitting to the broadband photometry, still makes this 
source a very luminous quasar with Mi=— 31.0 and changes 
our conclusions little. In Banerji et al. (2012, in preparation) 
we present follow-up observations of this source at multiple 
wavelengths in order to better characterise its nature. 

In the right panel of Fig. IIOI we plot the extinc- 
tion corrected absolute magnitudes versus redshift for all 
our reddened quasars co mpared to th e SDSS optically se- 
lected population of Shen et al.l (|201ll). The total reddened 
quasar sample from Paper I and this work covers an area of 
955 deg^ . We also list these extinction corrected i-band ab- 
solute m agnitudes in Table [6[ We have used a correction of 
+0.596 (JRichards et al.1 120061 ) to convert between Mi[z=2] 
and Mi for the SDSS quasars. The plot illustrates both the 
reddened absolute magnitudes as well as the de-reddened 
absolute magnitudes of our sample in order to demonstrate 
the effect of extinction on the inferred optical brightness. As 
can be seen from the figure, the extincted optical magnitudes 
lie well below the SDSS flux limit for this sample but cor- 
recting for this dust extinction makes the reddened quasars 



as luminous as some of the brightest SDSS sources and in 
some cases, brighter than any known SDSS quasar. This is 
to be expected given that we are selecting red quasars in a 
NIR flux-limited sample. Consequently, for a constant NIR 
flux, more heavily obscured sources are going to be more 
intrinsically luminous. This can be seen in the right-panel of 
Fig. [To] which also shows our i^-band selection boundary in 
the Mi — z plane. Our {J — K) > 2.5 colour selection broadly 
corresponds to selecting quasars at z^2 with E{B — V) > 0.5 
(Paper I), albeit with some scatter. We note however that 
while the NIR selection allows us to detect the most heav- 
ily obscured intrinsically luminous quasars at z ^2, it is 
biased against heavily obscured quasars with more mod- 
est luminosities at the same redshifts - i.e. a quasar with 
E{B — V) = 1.5 for example, would have to be one of the 
most luminous quasars known to make it into our K-band 
sample at z ~2. We discuss this selection effect further later 
in this Section. 

Fig. [To] clearly demonstrates that the reddened quasars 
we have found are some of the most intrinsically luminous 
quasars known at z^2. How does their space density com- 
pare to the optically selected population? Estimates of the 
space density of optical quasars in the SDSS for example 
are complicated by the fact that the SDSS is highly incom- 
plete in the redshift range 2.2< z <2.7, due to the similar- 
ity of Galactic star and qu asar colours at these redshifts. 
[Richards et al.l (|2002l , l2006l l have however studied this in- 
completeness in detail and corrected for it when deriving 
the quasar luminosity function. The selection functions de- 
rived by these authors were based on comparisons to sim- 
ulations but have recently been verifi ed empirically using a 
large sample of KX-selected quasars (|Maddox et al.ll2012l V 
We can therefore use the SDSS luminosity function to pre- 
dict the number of UV-bright quasars expected in the red- 
shift range 1.2<z<1.75 or 2.0<z<2.7 over our survey area of 
955 deg^ . These numbers of UV-bright quasars should then 
be directly comparable to the numbers of NIR bright red- 
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Figure 10. Left: Observed VLT-SINFONI spectrum for ULASJ1234+0907 (black) compared to the best-fit reddened model SED (red) 
with Ay=6.05 and the unreddened SED (blue) which has been corrected for extinction using the adopted reddening law fSection l3.1|) . We 
also plot a model with Ay =4.96 as the red dashed line. This model produces a slightly better match to the spectrum continuum than the 
A\/=6.05 model fit to the photometry. The UKIDSS H and X-band fiuxes that have been used to derive the higher Ay value are shown 
as the filled circles. The optical SDSS i-band bandpass as well as the near infra-red ULAS X-band bandpass are also shown for reference. 
ULASJ1234+0907 is a very extreme object with an effective reddening of 17.5 mags in the SDSS i-band (~2250A in the rest-frame) 
making it one of the most intrinsically luminous quasars known. Right: Absolute i-band magnitude versus redshift for our reddened 
quasars compared to the optical quasar population from SDSS DR7. The open circles correspond to the reddened absolute magnitudes 
and the filled circles represent the extinction corrected absolute magnitudes with the red dashed lines indicating the effect of extinction 
for each object. On correcting for the large extinction values, the reddened quasars mostly lie at the most luminous end of the optical 
sample with no spectroscopically confirmed SDSS quasars that are intrinsically as bright as ULASJ1234+0907 and ULASJ1539+0557. 

The plot also shows the tracks of a X = 17 quasar with different values of E{B — V) in the Mi z plane corresponding to the fiux-limit 

of our survey. 



Table 6. De-reddened absolute i-band magnitudes for all our 
spectroscopically confirmed red quasars. The quasars marked 
with an asterisk extend to a slightly fainter fiux-limit of 
Kas <18.9. They were targetted due to the presence of ancil- 
lary mult i- wavelength information but do not represent a com- 
plete sample and are therefore excluded from the space-density 
calculations. 



Name 


Dereddened M^ 


ULASJ0016-0038 


-27.49 


ULASJ0041-0021 


-28.94 


*ULASJ0141+0101 


-27.78 


ULASJ0144+0036 


-28.26 


ULASJ0144-0114 


-28.85 


ULASJ0221-0019 


-27.52 


ULASJ1539+0557 


-31.08 


*ULASJ1002+0137 


-27.29 


ULASJ1234+0907 


-31.81 


ULASJ1455+1230 


-28.29 


ULASJ2200+0056 


-29.16 


ULASJ2224-0015 


-27.96 



dened quasars in our sample. We note here that our sample 
extends to brighter absolute magnitudes than the SDSS but 
the parametric form of the luminosity function can reason- 
ably be extrapolated to these bright magnitudes in order to 
provide an estimate for the number of very luminous UV- 
bright quasars. The completeness corrected total number 
of UV-bright quasars at 1.2<z<1.75 or 2.0<z<2.7 over 955 
deg^ of sky predicted from the SDSS luminosity function, is 
given in Table O 

Our sample of reddened quasars also suffers from in- 
completeness. The first source of incompleteness arises be- 



cause not all the red quasar candidates in our UKIDSS 
sample were targetted for spectroscopy. There are six can- 
didates with colours consistent with identifications as red- 
dened quasars but without spectra. Best-fit E[B — V) values 
for these sources are calculated, just as for the spectroscop- 
ically confirmed sample, but assuming redshifts of 1.5 and 
2.5. Absolute i-band magnitudes can then be calculated as 
before. The results of this exercise are summarised in Table 
[71 It can be seen that the absolute i-band magnitudes are 
similar regardless of the redshift. However, the red [H — K) 
colours of these sources strongly favour the higher redshift 
assumption. 

The second source of incompleteness is due to morpho- 
logical misclassification of candidates, whereby some red- 
dened quasars, classified as galaxies in the i^-band, are not 
included in our sample. In Paper I we have shown that at 
the brighter end of our survey the maximum morphological 
incompleteness in our redshift range is ^20 per cent. 

Finally, as already pointed out in Fig. [TOl our flux- 
limited NIR sample only allows selection of the most intrin- 
sically luminous quasars ai z ^2 with large amounts of dust 
extinction. More modest luminosity quasars with similar or 
higher levels of dust obscuration require a deeper i^-band 
survey to be found so we could potentially be missing large 
numbers of highly reddened quasars at modest luminosities. 
We estimate our survey incompleteness to a typical Aa/=2.5 
quasar by assuming that every SDSS quasar at2.0<^<2.5 
for example, has this level of dust extinction, and counting 
the fraction of these dusty quasars that satisfy our survey 
selection criteria {K < 16.5 and {J — K) > 2.5; Vega) in 
each of the absolute magnitude bins in Table (8] We find 
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Table 7. Six quasar candidates for which no spectra were obtained along with their best-fit E{B - 
as well as de-reddened absolute i-band magnitudes corresponding to these two redshifts. 
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V) values assuming z =1.5 or z =2.5 



Name 


E(B-V),z=1.5 


Mi, z=1.5 


E(B-V), z=2.5 


Mi, z=2.5 


ULASJ0946+0016 


1.38 


-28.82 


0.80 


-28.59 


ULASJ1200+0423 


1.65 


-29.29 


0.98 


-29.05 


ULASJ1320+0948 


1.22 


-28.36 


0.69 


-28.12 


ULASJ0121+0107 


1.06 


-28.22 


0.60 


-28.04 


ULASJ0204+0106 


1.18 


-28.34 


0.67 


-28.13 


ULASJ2346-0038 


1.31 


-29.62 


0.75 


-29.38 



that our survey would recover >99% of Ay =2. 5 quasars 
with — 32 < Mi < — 30 at these redshifts with the complete- 
ness dropping to ~90%, 60% and 10% in the subsequent 
luminosity bins. We correct the numbers of quasars in our 
survey for all these sources of incompleteness to derive the 
upper limits quoted in Table (8] These are then compared to 
the numbers of UV-bright quasars from the SDSS luminos- 
ity function over the same redshift range and area, to derive 
an obscured fraction in each absolute magnitude bin. 



4.4 Is the Obscured Fraction Luminosity 
Dependent? 

Taken at face- value, Table [8] seems to suggest that the ob- 
scured fraction of broad-line quasars declines with intrinsic 
luminosity with these red quasars only making up a rela- 
tively small fraction of modest luminosity Type 1 quasars 
at z ^2. At very high luminosities however, heavily red- 
dened optically faint quasars dominate the population al- 
though with only small numbers at present to exhibit this 
trend. Is this trend physical or can it simply be explained 
by selection effects? While it is true that our survey is too 
bright to detect heavily reddened quasars with Ay > 3 at 
Mi > —28, should large populations of such sources exist, 
there are several reasons to suggest that this is probably not 
the case. Flat-spectrum radio-sour ce samples with very high 
quasar redshift completeness fe.g. ljackson et al.l ((2003)) in- 
dicate that the quasar population as a whole c annot consist 
of a very large fraction of obscured sources (JEllison et al.l 
[200l|). Although our sample size is still small, we have care- 
fully taken into account the various selection effects in our 
survey and we probe down to i^-band depths sufficient to de- 
tect most Ay =2. 5 quasars out to M^ < —28 at z^2. We con- 
clude that highly reddened quasars are likely to only domi- 
nate the counts of Type 1 quasars at the brightest absolute 
magnitudes, where NIR surveys may be more complete than 
optical surveys. 

The above arguments while still tentative due to the 
small numbers of quasars in our sample seem plausible based 
on physical arguments. The duty cycle for the obscured 
phase is likely to be longer for the most luminous quasars, 
consistent with a longer time required to assemble more mas- 
sive black holes. Seen another way, the obscured phase is 
more likely to be associated with the most intrinsically lumi- 
nous quasars with high mass accretion rates. These rapidly 
accreting, massive, obscured black-holes should preferen- 
tially lie at Z'-^2 corresponding to the main epoch of black- 
hole formation. Indeed, our two reddest and most luminous 
quasars ULASJ1234+0907 and ULASJ1539+0557 lie at the 
high redshift end of our sample at z ^2.5. Our selection 



criteria do not in any way discriminate against similarly lu- 
minous highly reddened quasars at 1.2< z <1.75 but both 
the quasars observed at these lower redshifts have compara- 
bly modest luminosities. Again, we caution against drawing 
strong conclusions from the very small numbers of sources 
in our sample but the fact that our most reddened and most 
luminous quasars are concentrated at the high redshift end 
of our sample, is also consistent with the most massive black 
holes assembling the bulk of their mass at early times via 
an obscured phase. 

Taken together, all these facts suggest that the obscured 
fraction may depend on quasar luminosity with a larger 
proportion of the most luminous quasars being preferen- 
tially observed in an obscured phase at z^2. Our sample 
is still however highly incomplete for obscured quasars at 
Mi > — 28 and a deeper K-hand sample is clearly needed to 
fully test this theory. 



4.5 Comparison to the 2MASS Sample 

Having demonstrated that our population of intrinsically lu- 
minous dust-reddened quasars at z^2 with E{B — V)> 0.5 
seem to represent a sizeable fraction of the total quasar 
population at very bright absolute magnitudes, we now 
make some explicit comparisons to previous studies of such 
objects using the sh allower (i^vega < 15) 2MASS survey 
(|Glikman et al.l 120071 ). Our survey represents a significant 
increase in depth compared to this previous work and ex- 
tends 1.5-2 mags fainter in the i^-band. Isolating only those 
quasars in the 2MASS survey that overlap the redshift range 
of our sample (1.2< z <1.75 or 2.0< z <2.7), there are 
only 13 reddened bro ad- line quasars over 2716 deg^ in the 
iGlikman et al.l (|2007l l sample. For each of these 13 QSOs, we 
use the K-band mag nitudes and dust extinctions given in 
iGlikman et al.l ((2003) along with the average {i — K) colour 
of our quasar models f^ l3.ip at that redshift, to calculate an 
absolute i-band magnitude. The counts are then scaled to 
our survey area of 955 deg^ and once again quoted in Table[8l 
Table [8] demonstrates that our NIR se lection uncover s much 
larger numbers of quasars at Z'-^2 than IGlikman et al.l (|2Q07l l 
but with similar intrinsic luminosities and much higher lev- 
els of dust-reddening. The z^2 2MASS sample has a median 
E{B — V^) = 0.4 while our sample over the same redshift 
range has a median E[B — V) — 0.8. 

We have also shown that the most heavily dust- 
reddened quasars in our sample lie at the highe s t reds hifts 
whereas the reddest quasars in IGlikman et al.l (|2007|) are 
preferentially at the low redshift end of their sample due 
to the shallow i^-band limit. If red quasars in which the 
red colours are attributed to dust extinction, represent a 
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Table 8. Number of quasars at 1.2<z<1.75 or 2.0<z<2.7 over 955 deg^ in bins of M^ for our survey of red quasars compared to the SDSS 
DR7 sample of blue quasars (corrected for completeness). The obscured fraction represents the fraction of reddened quasars relative to 
the total number of quasars in a particular absolute magnitude bin. In the case of our survey, the lower limits are computed from the 
raw counts of reddened quasars in this work while the upper limits take into account the various sources of incompleteness in the sample 
detailed in Section. 14.31 



M, 




N/955deg2(SDSS DR7) 


N/955deg2 (This Work) 


Obscured Fraction 


N/955deg2(F2M (Glikman07)) 


-32 to 


-30 


0.36 


2-3 


80-90% 


1 


-30 to 


-29 


3.3 


1-4 


20-60% 


0.7 


-29 to 


-28 


33 


4-17 


10-30% 


1.8 


-28 to 


-27 


304 


3-50 


1-10% 


1 



young phase in quasar evolution, they should preferentially 
be seen at z^2 corresponding to the peak of star formation 
and black-hole accretion. This is what is observed in our 
sample and suggests that other factors apart from dust ex- 
tinction in the quasar host may be r esponsible for the red 
colours in the reddest quasars in the iGlikman et alj (|2007l l 
sam ple. Selecting only qu asars with E{B — V) >0.5 from 
the IClikman et al.l (|2007h sample in order to overlap the 
typical reddening in the quasars in our sample, we find a 
median redshift of 0.6 for these quasars. At these redshifts, 
the host galaxy contribution to the quasar can be significant, 
affecting the quasar colours. In addition, at these lower red- 
shifts, the near infra-red J and iT-bands are also tracing 
hot dust emission. In a forthcoming paper, we demonstrate 
using WISE photometry for our quasars that the tempera- 
ture of this hot dust can significantly affect the colours of 
quasars at rest-frame wavelengths of 1-2/im independent of 
the amount of dust extinction i n the quasar host. The advan- 
tage of our sample over that of lClikman et al.l (|2007h is that 
at the higher redshifts, the hot dust emission is redshifted 
out of the NIR bands employed for the colour selection so the 
red colours can more safely be attributed to dust extinction 
in the quasar host. 



5 FUTURE WORK 

The current work is based on a sample selected using ^1000 
deg^ of data from the UKIDSS-LAS. The latest data release 
of the UKIDSS-LAS spans -3600 deg^ In addition, the on- 
going VISTA Hemisphere Survey will provide an additional 
~ 10,000 deg^ of high- latitude sky therefore easily allowing 
our sample sizes to be increased by an order of magnitude 
over the next five years as well as probing fainter fluxes in 
the i^-band. The combination of these near infra-red sur- 
veys with mid-infrared data from the all-sky WISE sur- 
vey will enable statistically significant samples of reddened 
Type 1 AGN to be assembled. We have already embarked 
on spectroscopic follow-up of these reddened AGN candi- 
dates selected from the initial data in the new surveys using 
Gemini-GNIRS and this new data will allow us to tighten 
our constraints on the exact obscured fraction of AGN at 
the main epoch of galaxy assembly and black-hole accretion 
as well as to probe in detail the dependence of this fraction 
on the quasar luminosity. 



6 CONCLUSIONS 

We have extended our search for NIR selected samples of 
reddened quasars to cover an area of almost 1000 deg^. 
Within this area, we have identified a total of 27 very bright 
K-hand sources that also have unusually red colours. Ten 
of these have been presented in Hawthorn et al. In this pa- 
per we present spectroscopic observations of another 13 can- 
didates. Five of these are confirmed to be highly reddened 
Type 1 AGN increasing our total reddened NIR bright AGN 
sample to 12. 

We combine the new sample and the sample from 
Hawthorn et al. in order to examine the properties of the 
optically obscured luminous Type 1 AGN population. Un- 
like previous studies, our sample relies on a homogenous 
selection at NIR wavelengths without needing detections in 
the optical and/or the radio. We carefully account for the 
effect of Ha equivalent widths on the infra-red colours of our 
quasars and by probing luminous sources at z >1, eliminate 
any contribution to their colours from the host galaxy. In 
particular, we reach the following conclusions: 

• Our sample of spectroscopically confirmed reddened 
quasars with redshifts all possess broad Ha emission indi- 
cating that they are canonical Type 1 AGN. 

• Almost all Ha line profiles are best fit by a broad 
(>5000kms~^) and intermediate (>1500kms~^) component 
and some of the lines show significant velocity offsets be- 
tween the two components indicating the presence of strong 
outflows in the broad-line region of the AGN. Such outflows 
are not commonly seen to affect the Ka lines in UV-luminous 
quasars and provide evidence that the reddened quasars are 
observed in a phase when they are in the process of expelling 
their gas and dust. 

• We use the broadband colours to fit model SEDs to all 
quasars with redshifts and find that significant dust extinc- 
tion of Ay =2-6 mags is required to fit the observational 
data making these quasars highly obscured at optical wave- 
lengths. The dust extinctions in our sample are comparable 
to those seen in the submillimetre galaxy population at sim- 
ilar redshifts and this level of dust obscuration pushes the 
quasars well below the flux limits of large-area optical sur- 
veys like the SDSS. 

• We use the line profiles to derive virial black-hole 
masses and Eddington fractions for all our reddened quasars. 
The mean virial black-hole mass is 2xlO^M0 and the mean 
Eddington fraction is 0.5. The reddened quasars are there- 
fore already as massive as the most massive optically se- 
lected quasars but still accreting mass at a relatively high 
rate. We argue that they probably represent a short-lived 
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phase in quasar evolution when the most massive starbursts 
are transitioning to UV-luminous quasars. 

• Two of our confirmed reddened quasars, 
ULASJ1539+0557 {z = 2.658) and ULASJ1234+0907 
(z = 2.503) are intrinsicahy brighter than any known SDSS 
spectroscopic quasars and have mass accretion rates in 
excess of lOOMoyr"^. They both he at the high-redshift 
end of our sample consistent with a picture where the most 
heavily reddened luminous quasars are accreting the bulk 
of their mass at early times via an obscured phase. 

• We consider the possibility of lensing in our population 
and conclude that while we cannot rule out the possibility 
that individual objects are lensed, our sample is generally 
fainter in the K-hand than all the known lensed reddened 
quasars previously found. Crucially, many of our sources are 
undetected in the SDSS i, z and UKIDSS Y and J bands 
where a lensing elliptical galaxy at z ^1 should be visible. 
The bulk of our reddened quasar population is therefore 
unlikely to be lensed. 

• The broad lines combined with the large dust extinc- 
tions observed in our sample strongly point to a specific 
phase in the quasar life cycle, rather than orientation, be- 
ing responsible for the red colours. We conclude that our 
quasars are likely to be observed in a phase of coeval galaxy 
and black-hole growth when the starburst is decaying but 
the dust has not yet fully cleared preventing these sources 
being selected as optical quasars. 

• We use our sample to place new constraints of the frac- 
tion of obscured quasars and find tentative evidence that 
this is a function of intrinsic luminosity and dust reddening. 
Although a large population of modest luminosity heavily 
dust-obscured quasars is not ruled out by our data given 
the K-hand flux limit of our sample, other independent lines 
of evidence point to this not being the case. At bright abso- 
lute magnitudes however, where we are reasonably complete, 
the obscured quasars seem to make up a sizeable proportion 
of the Type 1 quasar population suggesting that NIR sur- 
veys may be more complete than optical surveys at the very 
bright-end of the quasar luminosity function at z ^ 2. 

• The large obscured fractions in the most intrinsically 
luminous quasars are consistent with these sources having a 
longer obscured duty cycle. Seen another way, the obscured 
phase of a quasar is probably likely to be associated with 
its most intrinsically luminous stage as it corresponds to a 
time when the black-hole is accreting mass very quickly. 

• Given these results and their possible implications for 
quasar evolution, a larger sample of K-hand selected quasars 
pushing to fainter K-hand magnitudes is clearly required 
to test the hypothesis put forward here that the obscured 
fraction of Type 1 AGN is a function of quasar luminosity. 

In conclusion, we have presented a NIR-selected sam- 
ple of reddened quasars that have been chosen to lie at 
1.5< z <3 corresponding to the peak of star formation and 
black-hole accretion in the Universe. Unlike previous stud- 
ies, we detect populations of quasars that lie well below the 
SDSS detection limits as well as radio-quiet quasars that are 
not present in the FIRST radio survey, thus demonstrating 
the potential of new large-area infra-red surveys in studying 
the dust-obscured quasar population. We find that these lu- 
minous dusty quasars are likely being observed in a phase 
of coeval galaxy and black-hole growth and are in the pro- 



cess of expelling their gas and dust before turning into UV- 
luminous quasars. Our ongoing follow-up of these sources at 
long wavelengths will allow their host galaxy properties to 
be constrained for the first time. 

Having demonstrated the success of NIR surveys in se- 
lecting reddened quasars, we have embarked on a spectro- 
scopic campaign to follow-up new sources being discovered 
in the VISTA Hemisphere Survey and the Wide Infrared 
Survey Explorer. These new large area datasets will enable 
much larger samples of dust obscured quasars to be assem- 
bled over the coming years and push to fainter K-hand lim- 
its thereby allowing accurate constraints to be placed on 
the physical properties and global space densities of heavily 
obscured broad-line quasars. 
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